A Search for Sub-TeV Gamma Rays in Coincidence with Gamma Ray Bursts 
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We report on a study of sub-TeV {E-, > 10 GeV) gamma-ray-induced muon secondaries in co- 
incidence with BATSE gamma ray bursts (GRBs). Each TeV gamma ray striking the atmosphere 
produces ~0.2 muons whose identity and angle can be measured by the Project GRAND array. Eight 
GRB candidates were studied; seven were selected based upon the optimum product of (detected 
BATSE fluence) x (GRAND's acceptance). One candidate was added because it was reported as 
a possible detection by the Milagrito collaboration. Seven candidates show a positive, though not 
statistically significant, muon excess. The only significant possible coincidence shows an excess of 
466 ± 171 muons during the BATSE T90 time interval for GRB 971110. The chance probability of 
such an excess in GRAND's background at the time of this event is 3 x 10~^. The chance probability 
of observing such an excess in one of the eight bursts studied here is 0.025. If this event is real, the 
implied fluence of energetic (> 10 GeV) gamma rays necessary to account for the observed muon 
excess would require that most of the GRB fluence arrived in the form of energetic gamma rays. 

PACS numbers: PACS Numbers: 95.85.Ry, 98.70.Rz, 98.70.Sa, 14.60.Ef, 95.55.Vj 



I. INTRODUCTION 

The mystery of the astrophysical origin for gamma ray 
bursts (GRBs) has been with us for some time. As of yet 
there is no consensus explanation for them. Neverthe- 
less, a likely scenario is a burst environment involving 
collisions of an ultra relativistic e"*" — e~ plasma fire- 

ball 0, i a. These fireballs may produce low-energy 

gamma rays either by "internal" collisions of multiple 
shocks P, a , or by "external" collisions of a single shock 
with ambient clumps of interstellar material 0. 

In either of these possible paradigms, however, it seems 
likely that energetic ('^TeV) gamma rays and/or neu- 
trinos might also be produced along with the low en- 
ergy GRB. For example, inverse Compton scattering of 
ambient photons with relativistic electrons could pro- 
duce high-energy gamma rays 01 ■ Alternatively, baryons 
would be accelerated along with the pair plasma to very 
high energies ^, 9 , lOj . Synchrotron emission from en- 
ergetic protons 0jOr hadroproduction of pions in the 
burst environment [lO| and subsequent tt" gamma ray de- 
cay might also yield a spectrum of energetic gamma rays. 

Thus, it is plausible that energetic gamma rays could 
be emitted in coincidence with a lower-energy GRB. 
These energetic gamma rays would, however, be attenu- 
ated by the infrared background. Hence, their detection 
will require that the source be nearby and/or that their 
flux constitute a significant fraction of the energetic out- 
put of the source. Nevertheless, such energetic gamma 
rays, if detected, could provide valuable clues as to the 
baryon loading, Lorentz factor, and ambient magnetic 
field of the relativistic fireball. They might also provide 
a means to distinguish between an internal vs. external 
shock origin for the bursts. This paper, therefore, reports 
on an independent search for the possible coincidence of 
high-energy gamma rays with GRBs. 

There presently exists at least some evidence for a pos- 
sible association of energetic gamma rays with low-energy 



GRBs in previous literature. EGRET detected seven 
GRBs which emitted high energy photons in the ^100 
MeV to 18 GeV range [nUIllia. There have also been 
some results from the Tibet air shower array suggestive 
of gamma rays beyond the TeV range 0, ^3 j although 
these results were not claimed as a firm detection. There 
has also been reported evidence for TeV emission in one 
burst out of 54 BATSE GRBs in the field of view of 
the Milagrito detector |^^. The chance probability for 
that event (out of 54 trials) was found to be 1.5 x 10~^. 
During the present study, this Milagrito event was also 
within the field of view of GRAND, but at a relatively 
low elevation. 

In this paper we report a marginal {2.7a) detection 
of energetic (E^ > 10 GeV) gamma rays in coincidence 
with one BATSE GRB (971110). The chance probability 
of such an excess from this burst is 3 x 10~^. The chance 
probability of such an excess from our sample of eight 
bursts is 0.025. The constraint that this detection might 
place on the spectrum of energetic photons from the burst 
environment is discussed. The analysis of GRAND's data 
for the GRB reported as a possible detection by Milagrito 
finds an insignificant excess of muons (1.2(T, including 
statistical and systematic errors). 



II. PROJECT GRAND 

GRAND is located just north of the University of Notre 
Dame campus, approximately 150 km east of Chicago 
and 220 m above sea level at 86.2° W and 41.7° N. It de- 
tects cosmic ray secondaries at ground level by means 
of 64 tracking stations of proportional wire chambers 
(PWCs)[l3. Each station has four PWC detectors; each 
detector contains an x-plane (x=eastward) and a y-plane 
(y=northward) oriented horizontally and stacked verti- 
cally (z) The planes have an active area of 1.29 
m^ comprised of 80 detection cells (each 14 mm x 19 mm 
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X 1.1m). Each secondary muon is measured to 0.26° ab- 
solute precision (average value in each of two orthogonal 
planes) |23|. At present, GRAND 's total detector area 
is 83 m^; earlier data had smaller areas. A 51 mm thick 
steel plate inserted between the third and fourth PWCs 
allows muon tracks to be distinguished from electrons; 
96% of muon tracks are correctly identified as muon 
tracks. Since, for single tracks, only 1/4 are electrons 
and 4% of these electrons are misidentified as muons, the 
muon sample has only 1% contamination from electrons. 
The data presented here are from single-track triggers 
which ignore time coincidences between stations; only 
muon candidates are selected. 

GRAND utilizes the fact that gamma rays have a 
detectable signal of muons from gamma-hadro produc- 
tion in the atmosphere. The pions thus produced subse- 
quently decay to muons which can reach detection level 
making it possible to study coincidences between GRBs 
and gamma ray showers in the Ey > 10 GeV energy re- 
gion. This energy threshold (10 GeV) depends slightly 
upon the spectral index of the gamma ray spectrum (see 
Figure 6 in Ref. [23) and is not a sharp threshold. It 
has been estimated that the GRB rate for a thresh- 
old energy larger than 200 GeV is ~10 GRBs per year; 
for the ~ 10 GeV threshold energy of GRAND, this rate 
should be even greater. 

Since the study reported here is similar to and in- 
cludes the event reported by the Milagrito collaboration, 
we note the differences between GRAND and Milagrito. 
Whereas GRAND detects secondary air-shower muons 
with PWCs, Milagrito detects secondary air shower 
charged particles by Cerenkov light as the shower par- 
ticles traverse a light tight water reservoir. The main 
differences between the two detectors are that Milagrito 
has a larger active detection somewhat better an- 

gular resolution in the single-track mode of interest here 
(-- 1° for Milagrito vs. ~ 5° for GRAND), and a higher 
detection threshold (~ 1 TeV for Milagrito vs. ^ 10 GeV 
for GRAND). The lower detection threshold for GRAND 
implies that it is sensitive to a lower energy part of the 
primary gamma ray spectrum which is not as likely to 
have been extinguished by internal or intergalactic 
absorption [2^ 01 • Thus, even though no individual 
detection was overwhelmingly significant, the fact that a 
slightly positive signal to background was deduced for all 
but one of the eight candidate bursts investigated might 
suggest that >10 GeV emission in association with low- 
energy GRBs is not altogether uncommon. 

GRAND is a unique detector facility which measures 
the angles of single tracks and identifies which are muons. 
Inferring the implied gamma ray flux from the detected 
muons, however, requires confidence in the ability to sim- 
ulate the muon production from gamma hadroproduction 
and muon propagation in the atmosphere. Recently, we 
have made a detailed analysis [2ll l25l | of the spatial and 
energy distribution of muons in 7-induced air showers. 
These simulations are based upon the FLUKA Monte- 
Carlo (MC) code. Unlike most MC codes used in cos- 



mic ray research, this simulation is not specialized for 
this particular field but is a multipurpose particle trans- 
port code which has been tested in many diverse appli- 
cations such as proton and electron accelerator shield- 
ing, calorimetry, medical physics, etc. It has therefore 
been verified against a large amount of nuclear experi- 
mental data and indirectly validated by comparisons with 
shower measurements obtained both at accelerators and 
in cosmic ray experiments. Of particular relevance to the 
present study is the fact that this code has been shown 
|26j to accurately predict hadron-generated muon spectra 
at different heights in the atmosphere. Hence, we expect 
its application here to photon-generated muon spectra be 
good to the ^few % statistical accuracy of the simula- 
tions. 

The MC simulation of (25|] shows that a 1 TeV gamma 
ray normally incident upon the earth's atmosphere pro- 
duces an average of 0.23 muons which reach GRAND. 
GRAND thus paradoxically uses muons as a signal for 
gamma ray primaries. In the energy region > 10 GeV, 
the muon statistics are quite high; the current all-sky 
rate for recording identified muons is about 2400 Hz or 
8.6 million muons per hour. 

These secondary muons are primarily the result of in- 
teractions of primary cosmic rays with the atmosphere 
producing pions, which then decay to muons. The pi- 
ons are produced at small angles relative to the primary 
cosmic rays; the non-interacting pions decay to muons 
which emerge at small angles relative to the pion; the 
muons are then deflected in the earth's magnetic field 
and scattered in the atmosphere resulting in an effective 
net angular resolution of about ± 5° (for the primary 
cosmic ray in each of the two orthogonal directions; this 
resolution depends slightly upon the primary energies or 
spectral index). The muon threshold detection energy is 
0.1 GeV for vertical tracks to penetrate the 50 mm steel 
plate; however, these muons must have been born with at 
least several GeV of energy to penetrate the atmosphere 
in order to reach the detectors. 

The precise response to primary gamma rays is de- 
scribed with a Monte Carlo calculation (see Figure 6 of 
j27l|): the result depends upon the assumed spectral in- 
dex. The primary gamma-ray spectrum is assumed to be 
of the form E^^ with a spectral index /3 = —2.41 (the av- 
erage of the spectral indices reported in the third EGRET 
catalog 13). This spectrum of primary gamma rays is 
then multiplied by the number of muons per gamma ray 
which reach detection level |2^ |23| . This determines the 
number of detectable muons as a function of the primary 
gamma-ray energy. Qualitatively, each primary gamma 
ray produces a number of muons. This number at first 
increases sharply for gamma-ray primary energies from 
1.5 to 10 GeV and then falls slowly for energies above 
10 GeV. For a softer spectral index, the muon response 
peaks at a slightly lower primary gamma-ray energy and 
then falls off more rapidly above 10 GeV. For conve- 
nience, we characterise this response shape as having a 
threshold of ^10 GeV for the primary gamma-ray spec- 
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tral indices of interest. 

grand's ability to correlate short bursts of muons 
with an identifiable source of primary radiation has been 
shown in a detection which was coincident with a solar 
flare on 15 April 2001. The statistical significance of this 
observation was at the level of 6 u for a ground level event 
of 0.6 hours duration [2^ . 



III. DATA ANALYSIS 

To analyze GRB coincidences the complete GRB ta- 
ble, the flux table, and the duration table were down- 
loaded from the BATSE archive 0. Two of the can- 
didate GRBs are listed in the BATSE 4b catalog, the 
other more recent events are in the BATSE archives. Ta- 
ble 1 lists some of the BATSE data including: the date of 
the trigger (GRB), the trigger number (Trig); the time 
duration for 90% of the burst's counts to occur (T90) 
in seconds, the right ascension (RA), declination (Dec), 
and the BATSE angular error (56), all in degrees. Next 
to these BATSE entries are the calculated angular eleva- 
tion (in degrees) above GRAND's horizon (denoted Elev) 
and our selection criterion (LogLk) described below. The 
last three columns of Table 1 summarize the muon secon- 
daries observed by GRAND within a ±5° square angular 
window centered on the burst location during the BATSE 
T90 time interval. 



A. Selection Criterion 

For each GRB, an approximate total acceptance factor, 
A, was calculated for the detector stations of GRAND: 

A~ [(1 — 5 X tan0^)(l — gx ta,ndy) x cosfl] x cos^ 6 (1) 

where the geometrical factor g = d/L, with d the ver- 
tical spacing between the top and bottom plane (0.61 
m), and L the length of a single plane (1.1 m). The an- 
gles 9x and 9y are the angle of the track from vertical 
{6, the complement of Elev) projected onto the xz- and 
yz-planes, respectively. The quantity in square brackets 
is the geometrical acceptance of a detector station; this 
acceptance has been multiplied by cos^ 9 to account for 
the added absorption of a muon traversing an increased 
path length through the atmosphere for tracks inclined 
from the vertical. 

As a rudimentary criterion to select which events to 
analyze, the relative likelihood LogLk of GRAND to ob- 
serve each GRB, is taken as proportional to the (base 10) 
logarithm of the product of the total acceptance factor, 
equation (1), times the BATSE fluence in the highest of 
four energy bins (i.e. > 300 KeV). The GRBs with 
the 21 largest values for LogLk > 3 from the BATSE 
archive were selected for further consideration. In ad- 
dition, it was recognized that the Milagrito event was 
in the field of view of GRAND, and even though it did 



not satisfy the selection criterion, it was included in the 
analyses. 

The entries in Table 1 are ordered with the highest 
value of LogLk at the top. This LogLk factor, how- 
ever, did not take into account the fact that the array 
was under continuous construction during this time and 
had varying (usually increasing) numbers of operational 
detectors at a given time. The data for each GRB an- 
alyzed were checked to ensure no huts were turning on 
or off during the time of analysis for that GRB. Data 
for 11 of the top 21 GRBs by this LogLk criterion were 
available on archived data tapes. Of the 11, one was 
found to have individual stations with large time depen- 
dent inefRciencies and was discarded. In addition, three 
of the candidate GRBs were before 1994 when the detec- 
tor area was small and had errors in the clock's seconds 
bit. Since these three bursts had the shortest T90 times, 
they were most sensitive to the precise time. Further- 
more, the BATSE angular errors were largest (two were 
comparable to the total width of our analysis window). 
Therefore, these three events were eliminated from the 
present analyses. [These three events were included in a 
previous preliminary analysis [3H with chance probabil- 
ities of -f 1.6, -0.1, and -1.4cr (Stat.)] The calculation of 
time from the data was checked to ensure that the times 
of the remaining GRBs are correct. 

In a preceding preliminary analysis [3l| , we employed 
an angular window using fixed angles of right ascension 
(RA) and declination (Dec). However, the signal and 
background data in this analysis showed systematic vari- 
ations due to: 1) The local angles 6x and dy change 
with time during a burst, which means that event tracks 
are recorded with differing numbers of wire combinations 
leading to systematic fluctuations in the comit rate. 2) 
The detector's solid angle increases as a fixed RA and 
Dec moves toward the vertical (or conversely). 3) There 
is dead time (tens of milliseconds) whenever: a) the data 
input is stopped in order to write an event buffer to mag- 
netic tape; or b) a trigger is received for a concurrently 
running air-shower experiment. These systematic effects 
were accounted for in the current analysis method as fol- 
lows: 

First, the BATSE GRB locations in RA and Dec were 
transformed to a horizon coordinate system of elevation 
and azimuth and then projected onto the xz-plane {9x) 
and the yz-plane {9y). A window of ±5° in A9x and A9y 
was centered on the location of the GRB. 

As time progressed during the T90 burst interval, the 
local windows of 9^ and 9y stay fixed thus removing prob- 
lems (1) and (2) above. To correct for the dead time (3), 
the total event rate over the whole sky was employed as 
a high-statistics measure of the live time of each time 
bin; each bin's data were corrected for its corresponding 
live time. After this correction, the data inside the an- 
gular cuts were almost independent of the random dead 
times. These corrected numbers of muons inside the an- 
gular window were studied to obtain the (background -I- 
signal) within the T90 interval. The background was de- 
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termined during a time interval of «20xT90 before the 
start of the BATSE trigger (except in the case of GRB 
971110 it was wlOxTQO as the longer interval of time 
was not available on the same data tape for this event). 



B. Significance 

The signal (Sig) calculated for a GRB is the differ- 
ence between the total counts Non inside the T90 inter- 
val (corrected for dead time) and the background counts 
Nof f normalized to the live time of the T90 time interval 
and corrected for dead time. The statistical significance 
Sig/SSig (number of standard deviations above back- 
ground) of each event was determined according to the 
likelihood ratio method of Li and Ma l32l , 



Sig/SSig = \/2l^Nonl'r 



1 + A / iVo, 



+ Nof fin 



(1 + A) 



A \Non+Noff 

Noff 



Non + Nof f 



1/2 



(?) 



The only possibly significant single observation among 
the eight GRBs in Table 1 corresponds to GRB 971110. 
The time distribution of the data before and after the 
BATSE event trigger are illustrated in Figure^ Here the 
data are plotted in a manner similar to the way in which 
the Milagrito event was presented, i.e., with two dif- 
ferent values for the time-bin resolution and for differ- 
ent durations before and after the GRB. In this figure 
negative background-subtracted events are also included 
giving added information on the fluctuations in the back- 
ground. There are excess events in the T90 window (and 
possibly before and after as well). The muon excess in 
the T90 window is 3.3cr above the a/TV fluctuations ex- 
pected from pure counting statistics. Nevertheless, it is 
also possible that the background fluctuates in excess of 
the expected \/N statistics. This could happen due to 
real variations in the muon arrival rate at a given angle 
or to variations in the detector response not already ac- 
counted for in our dead time correction. In either case, 
the additional systematic error associated with such pos- 
sible excess fluctuations can be deduced directly from the 
observed data rate as we now describe. 



where A is the ratio A = ton /toff- This method makes a 
best accounting of the true significance of an event when 
the background and event time intervals differ, i.e. A 7^ 1. 
For aU of these events A « 0.05 except for GRB 971110 for 
which A ~ 0.10. The statistical errors from this analysis 
are listed in the last column of Table 1. 

The data for these GRBs were tested before, during, 
and after the signal region to see if any station's data had 
excessive noise or erratic time dependences which could 
erroneously mimic a GRB signal. These rates exhibited 
no pathological time structure. 
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FIG. 1: Background-subtraced event rate in the Project 
GRAND array before and after the BATSE 6472 event trigger. 
The T90 interval for this burst is indicated by a horizontal line 
above each of the two histograms. 



C. Systematic Error 

As a check on systematic errors in the signal for GRB 
971110, similar angular sections of the sky (which have 
the same absolute values of 9x and 9y and thus the same 
average counting rates) and the same time interval (T90) 
but at different, neighboring times were analyzed. These 
time intervals span ±26 hours relative to the BATSE 
trigger for GRB 971110. The data on the tape were also 
analyzed for T90 intervals beginning at times before and 
after the GRB, and for different 9x and 9y locations in 
the sky but at the same elevation angle. This analysis 
was done in the same way as described in the preceding 
paragraphs. 

Figure |21 shows the distribution of fluctuations above 
and below background for 1587 separate analyses (not 
including the time and angle of the GRB 97110 event). 
The resulting distribution is centered on zero, and is ap- 
proximately Gaussian in shape. However, this distribu- 
tion has a larger width than that calculated from the 
Poisson random statistics of the counts within the T90 
and background intervals. The standard deviation width 
of this distribution is 171 counts, whereas the expected 
statistical deviation is only 141 counts based upon the 
background count rate. A total statistical plus system- 
atic standard deviation is 171 counts for this distribution. 
Adopting this total error as the quadrature sum of the 
statistical and systematic errors, then 97 counts are as- 
cribed as the additional systematic error in our analysis 
for GRB 971110. With this additional systematic error, 
the ratio of signal-to-noise becomes Sig/SSig = 2.72. 

The probability of a +2.7a fluctuation in a Gaussian 
white-noise distribution is 4.9 x 10~^. However, the his- 
togram in Figure 121 can be used to measure the probabil- 
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FIG. 2: Distribution of fluctuations above and below back- 
ground for 1587 random T90-intervals observed by GRAND 
near the time of GRB 971110 

ity of grand's detection being real without assuming a 
distribution which is Gaussian. In this histogram there 
are 10 random fluctuations either > +466 or < -466. 
This corresponds to a probability of 3 x 10"'^ for a back- 
ground fluctuation to produce a > +466 signal. For one 
event out of the 8 analyzed to exhibit this much devia- 
tion would correspond to a random probability of 0.025. 
Thus, the statistics of this event are interesting but not 
compelling. 

The total statistical plus systematic error derived 
above is consistent with a white-noise fluctuation that 
would scale with the statistical error. In this case, it is 
reasonable to take the ratio of total to statistical error as 
a constant independent of the particular T90 interval in 
question. This constant factor (~ 1.2) can then be used 
to roughly estimate the systematic error for the remain- 
ing events. These remaining events (except one) indicate 
positive (though not significant) deviations. 

IV. DISCUSSION 

Accepting for the moment the proposition that the 
muon excess associated with GRB 971110 may be real, 
it is worthwhile to consider some of the implications of 
this event. First, we wish to consider the possibility of 
whether these events could be explained by a naive ex- 
tension of the low-energy burst spectrum. Indeed, it is 
known from EGRET detections uli that it is possible for 
the spectrum of the low-energy burst detected by BATSE 
to extend up to ~ 10 GeV with no change in the spectral 
index. On the other hand, such an extension is not guar- 
anteed. For example, no sub MeV photons were observed 
at the time that EGRET observed an 18 GeV gamma ray. 

To explore the possibility that the observed muon ex- 
cess could be explained by a straightforward extension 
of the BATSE spectrum, the GRB 971110 spectrum was 
fit using the standard BATSE empirical spectral model 



d^iV _ / ^ \" / E{2 + a) \ 
dEdAdt ^ ^ \im KeV ) ^"^^ \ Epeak ) 

for low energy, E < {a ~ [3)Epeakl{'^ + a), and 

\ [100 KeV(2 + a)] J ) ' 

for E > {a — l3)Epeak/{'2 + «)■ The best-fit parameters 
for this burst are a = -1.02 ± 0.04, (3 = -2.33 ± 0.11, 
B = 0.0095 ± 0.0003, and Epeak = 303 ± 17 KeV, with a 
reduced of 0.93. 

The solid line in Figure |21 shows the optimum fit to 
the BATSE spectrum extended to the energy range of 
GRAND. Note that this naive extension is probably an 
over estimate, since absorption due to pair production 
from interactions of gamma rays with the intergalactic 
infrared background is expected to become significant 
above ~ 200 GeV |23,|2|. 

This extended spectrum was folded with the known ef- 
ficiency, , of high-energy photon to muon conversion in 
the atmosphere (in the region of GRAND « 0.23E^-^y 
[2^) to yield the spectrum of muons at GRAND expected 
from this extrapolation of the BATSE spectrum (dotted 
curve in Figure|2Jl . The broad peak of the muon spectrum 
at around 10 GeV illustrates the peak (threshold) sensi- 
tivity of GRAND for this particular gamma ray spec- 
trum. Integrating this muon spectrum over the primary 
gamma ray energy, Ej , yields the number of muons per 
area per time based upon this extrapolation. Multiply- 
ing this result by the effective muon detection area of 
GRAND at the time of this burst and by its T90 time 
interval yields 0.3 ± 0.6 muons-well below the observed 
466 ± 171 excess muons detected for this event. Even 
if all of the fit variables are adjusted to their respective 
upper la limits, only four muons would be expected. 

Clearly, a simple, naive (very long) extrapolation from 
the BATSE spectrum is inconsistent with the observed 
muon excess at GRAND. In a related paper the spe- 
cific constraints which these data place upon possible 
mechanisms for energetic gamma ray emission such as en- 
ergetic proton-synchrotron emission inverse-Compton 
scattering from relativistic electrons, or hadronic produc- 
tion of pious in the burst [j^ are discussed. 



V. CONCLUSION 

We have completed a search for evidence for energetic 
sub-TeV gamma rays in coincidence with low-energy 
gamma ray bursts. No convincing evidence is found, 
though there is a possible 2.7a detection associated with 
the best candidate burst. There is an insignificant, 1.2a, 
muon excess observed in association with the Milagrito 
event (GRB 970417a), as well as a slight positive excess 
for all (but one) of the remaining events investigated. 
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FIG. 3: Extrapolation of the best-fit energy spectrum for 
BATSE GRB 971110. A portion of the BATSE data are 
shown in the inset. Multiplying the extrapolated BATSE 
gamma ray spectrum by the gamma to muon conversion effi- 
ciency (as calculated by FLUKA H^) gives the spectrum of 
muons which would be observed by GRAND. 

Based upon an analysis of the most significant event, 
we conclude that if the detected muon excess is real, then 



probably a new sub-TeV component to the gamma ray 
spectrum is required. A determination of the magnitude 
of this component, however, will require knowledge of the 
source spectrum as well as the effects of photon absorp- 
tion both in the burst environment and cnroute. Studies 
along this line imply that most of the GRB emission 
would be in sub-TeV gamma-rays. 
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TABLE I: Summary of Events Analyzed. 



GRB 


Trig 


T90 


RA 


Dec 


se 


Elev 


LogLk 


Nor, 




± aTot 


CStat 


971110 


6472 


195.2 


242 


50 


0.6 


81 


5.18 


18,286 


17,820 


466 ± 171 


141 




7343 


62.5 


229 


42 


0.4 


56 


5.13 


1 079 


1 076 


3 ± 36 


30 


940526 


2994 


48.6 


132 


34 


1.7 


66 


4.68 


498 


478 


20 ±28 


23 


980420 


6694 


39.9 


293 


27 


0.6 


68 


4.02 


1,456 


1,417 


39 ±47 


39 


960428 


5450 


172.2 


304 


35 


1.0 


70 


3.83 


3,990 


3.933 


57 ±78 


64 


980105 


6560 


36.8 


37 


52 


1.4 


79 


3.46 


2,214 


2,229 


-15 ±61 


50 


980301 


6619 


36.0 


148 


35 


1.3 


76 


3.17 


2,053 


2,007 


38 ±56 


46 


970417a 


6188 


7.9 


290 


54 


1.6 


62 


2.08 


186 


166 


20 ±17 


14 



NOTE: T90 is in seconds. Angles RA, Dec, 69, and Elev are in degrees. 



